REPUBLIQUE FRANCAIS 



PCT7EP200 5/ 0 1 302 T 

lm «Wlsd 



I INSTITUT 
NATIONAL OE 
LA PROPRIETE 
INDUSTRIELLE 



BREVET D' INVENTION 



CERTIFICAT D'UTILITE - CERTIFICAT D'ADDITION 

EPO-DG 1 

COPIE OFFICIELLE 0 7. 02. 2006 

; '(76) 



Le Directeur general de l' Institut national de la propriete 
industrielle certifie que le document ci-annexe est la copie 
certifiee conforme d'une demande de titre de propriete 
industrielle deposee a I'lnstitut. 

Fait h Paris, le L7.iDI.-JiS 



Pour le Directeur general de I'lnstitut 
national de la propriete industrielle 
Le Chef du Departement des brevets 




Martine PLANCHE 



INSTITUT 
NATIONAL D E 
LA PROPRIETE 
INDUSTRIELLE 



SIEGE 

26 bis. rue do Saint-Petcrsbourg 
75BOO PARIS cede* 08 
Tolophone : 33 (0)1 53 04 53 04 
Telscopie:33 (0)1 53 04 45 23 
wmv.inpl.fr 



ETABLISSEMENT PUBLIC NATIONAL CREE PAR LA LOI N» 51-444 DU 19 AVRIL 1951 




1er depot 



IMIIITI1ILU 



26 bis, we da Saint Petersburg • 75800 Paris Codex 08 
Pour vous informer : INPI DIRECT 

0,15€TTC/ma 

Telecople : 33 (0)1 53 04 52 65 



SSTWHOV 2004 
ueu 75 INPI PARIS 34 SP 

N« D'ENREGISTREMENT 04 1 2 4 7 1 
NATIONAL ATTRIBUE PAR L'INPI 

DATE DED£p6t ATTRIBUTE /j Q\J ^ ^QQ^ 

PAR L'INPI ^ 1 



Voft rof Frances pour ce dossier 

(facultalif) BFF040188/AH 



CERTIFICAT D'UTILITE 

Code da la proprleta Intollectuelle - Llvra VI 

REQUETE EN DELIVRANCE 
page 1/2 

Cet Imprlma est & remplir llslblernent a I'encre noire 



mm 



N° 11354-03 



BR1 



DB 540 OW/ 030103 



|] N0M ET ADRESSE DU DEMANDEUR 0U DU MANDATAIRE 
A QUI LA CORRESPONDANCE DOIT ETRE ADRESSEE 

■ 

CABINET PLASSERAUD 
65/67 rue de la Victoire 

75440 PARIS CEDEX 09 



Confirmation d'un depdt par telecopio 



3 NATURE DE LA DEMANDE 



Demande de brevet 



Demande de certificat d'utilite 



Demande dlvlslonnaire 

Demande de brevet initiate 
ou demande de certificat d'utilite' initiate 



Transformation d'une demande de □ 

brevet europeen Demande de brevet in itiale N° 

TITRE DE L' INVENTION (200 caracteres ou espacoi maximum) 



□ N° attribue par 1'INPl a la telecopfe 



Cochez I'uno des 4 cases sufvantes 



m 



□ 



Date 1 I I I I > I 1 
Date I I I I I ' I 1 



Date 1,1 1 I I I I > I 



Proc6de, installation et programme d'ordinateur pour estimer I'effectif initial d'une population d'acides nucleiques, 
notamment par PCR. 



DECLARATION DE PRIORITY 
OU REQUETE DU BENEFICE DE 
LA DATE DE d£p6T D'UNE 
DEMANDE ANT&RIEURE FRANQAISE 



Pays ou organisation 

Date I i 1 i I i i i I N° 
Pays ou organisation 

Date 1 i 1 i 1 i i i 1 N° 
Pays ou organisation 

Date I ■ I ■ I i i i 1 N° 

□ S'il y a d'autres prlorites, cochez la case ot utilisez I'iniprime wSujte» 



DEMANDEUR (Cochez Tune des 2 cases) 



Norn 

ou denomination sociale 



Prenoms 



Forme juridique 



N° SIREN 



1 L "1.J..I 



Code APE-NAF 



Domicile 
ou 



Rue 



Code postal et ville 



Pays 



Nationality 



N° de telephone (faculiatij) 



Adresse electronique (facultatif) 



Personne morale 



PI Pereonne physique 



B10-RAD PASTEUR 



J L 



3 Boulevard Raymond Poincare 



|9 i2 |4 1 3 )0 1 MARNES-LA-COQUETTE 



Francaise 



N° de telecopie (facultatij) 



[g] S'il y a plus d'un demandeur, cochez la case et utilisez I'lmprime «tSufte» 



RemnlirimperativementJa_e tt 



1er depot 




CERTIFICAT D'UTIIITE 

REQUETE EN DELIVRANCE 
page 2/2 




REMISE OES PIECES 

mte 24 NOV 2004 
UEU 75 INPI PARIS 34 SP 

N» O'ENREGISTREMENT Q4 12-471 

NATIONAL ATTRIBUE PAR L'INPI 



DBMOW/ 191203 



MANDATAIRE lieu) 



Nom 



Prenom 



Cabinet ou Society 



Nationality 



N °de pouvolr permanent et/ou 
de lien contractuel 



Adresse 



Rue 



Code postal et ville 



Pays 



N° de telephone (facultatij) 



N° de telecopie (facultatij) 



Adresse electronique (facultatij) 



E| INVENTEUR (S) 



CABINET PLASSERAUD 



65/67 rue de la Victoire 



17 i5 i4 4 iQ I PARIS CEDEX 09 



FRANCE 



01 40 16 70 00 



01 42 80 01 59 



Les inventeurs sont nficessairement des personnes physiques 



Les demandeurs et les inventeurs 
sont les mernes personnes 



□ Oui 

H Non ; 



Pans co cas rempHr le formulaire de Designation d'inventeur(s) 



RAPPORT DE RECHERCHE 



Unlquoment pour une demande do brevet (y compris division et transformation) 



£tablissement immediat 
ou etabllssement differe 



m 
□ 

Chobc a fair© obligatoirement au ddpot (cf. Notice explicative Rubrique S) 



REDUCTION DU TAUX 
DES REDEVANCES 



Uniquement pour tos personnes physiques 

□ Requise pour la premiere fois pour cette invention (joindre tin avis de non-imposition) 

□ Obtenue anterieurement a ce depdt pour cette invention (joindre une copie de la 
decision d'admission a ^assistance gratuite ou indiquer sa reference ) : AG | , , , , j 



Jg] SEQUENCES DE NUCLEOTIDES 
ET/OU D'ACIDES AMINES 



□ Cochez la case si la description contient une liste de sequences 



Le support electronique de donnees est joint 

La declaration de conformite de la tiste de 
sequences sur support papier avec le 
support electronique de donnees est jointe 



□ 
□ 



Si vous avez utilise I'imprimd «Suite», 
Indiquez le nombre de pages jointes 



1 PAGE 



SIGNATURE DU DESV1AWDEUR 

OU DU PJ3AMDATA1R£ 

(Nom et qualftd du signataire) 

Albert HASSINE 
N° 04-0603 




VISA DE LA PREFECTURE 
OU DE L'INPI 




La loi n°78-17 du 6 Janvier 1978 relative a Tinformatique, aux fichiers et aux liberies s'applique aux reponses faites a ce formulaire. 
Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de I'INPI. 




1er depot 




■ IK ft HUT 
NATIONAL CI 
LA rtOPBIUI 



jR6s8rv6a.'lNPI,. 



BREVET D'BNVENTION 

CERTIFICAT D'UTILIT^ N° U354-03 

Code do la proprlete Intellectuelle - Uvro VI 

REQUETE EN DELIVRANCE 

Page suite N° A. I A. 



/SUITE 



REMISJ 
DATE 



SSltoV 2004 
ueu 75 IN PI PARIS 34 SP 

0412471 



N* D'ENREGISTREMENT 
NATIONAL ATTRlBUt PAR L'INPI 



Cet imprimd est a remplir lisiblement d I'encre noire oa o w 7210103 



Vos r6f Srences pour ce dossier \facultalij) 


DCCH/tfH DO/A LJ 

or rU4U 1 00/ An 


O DECLARATION DE PRIORITY 
OU REQUITE DU BENEFICE DE 
LA DATE DE DltPfiT D'UNE 
DEMANDE ANTtRIEURE FRANQAISE 


Pays ou organisation 

Hate 1 1 1 1 1 1 1 1 1 N° 
rays ou organisation 

rtato It II 1 1 1 1 N° 

Pays ou organisation 

Date 1 1 1 1 1 1 1 1 1 N° 


Hf. DEMANDEUR (Cochoz I'uno des 2 cases) 


|X) Person ne morale [ J Perso nn a physique 


Nom 

ou denomination sociale 


CENTRE NATIONAL DE LA RECHERCHE SCIENTIFIQUE - CNRS - 


Prenoms 




Forme juridique 


Pt^Kliocomont Di ihlif* Qoiontifim ita *Tor > h\nnlonlnMO PPQT 
tZlauHSScinclH rUUIIO, OOlol lUUMUtJ ol 1 cLI H lUIULJIMUo CrOI 


N° SIREN 


1 1 1 1 1 1 1 I 1 I 


Code APE-NAF 


1 , I 1 1 


Domicile 
ou 

siege 


Rue 


*3 n io Mir*Hol Anno 
0, rue IVIILIIcl Milyc 


Code postal et viUe 


I7.5i7i9i4I PARIS Cedex 16 


Pays 




Nationality 


Frangaise .* • 


N° de telephone \facultaiifi 




N° de telecopie {facultatij) 


i " ! 


Adresse Electronique tfacultatjfl 




H DEMANDEUR (Cochez Tune des 2 cases) 


□ Porsonne morale □ Perconne physique 


Nom 

ou denomination sociale 




Prenoms 




Forme juridique 




N° SIREN 




Code APE-NAF 


IiiiI 


' Domicile 
ou 

siege 


Rue 




. Code postal et ville 


1 i i i i 1 


Pays 




Nationality 






N° de telephone \facultatij) 




N° de telecopie tfacuitatij] 




Adresse electronique {facultatij) 




E3 SIGNATURE DU DEMANDEUR A!bert HASS1NE A A {[ 

OU DU MANDATAIRE N° 04-0603 r Tpf 

(Nom et quality du olgnataire) M [J/^f* 

j 


VISA DE LA PREFECTURE 
JDU DE UIWPI 



l_<3 IUI II / VTA/ WU W J«I"'VI w . ~ ' ~ - - ' 

Elle garantit un droit d'acces et de rectification pour les donnees vous concernant aupres de 1'INPI 



1 er depot 



1 

Procede, installation et programme d 1 ordinateur pour 
es timer l'effectif initial d'une population d'acides 
nucleiques, notamment par PCR 

5 La presente invention concerne l 1 estimation de l'effectif 
initial d ? une population d ! int6ret dans un echantillon 
soumis a succession de reactions d 1 amplification . 

La pr§sente invention trouve une application avantageuse 
10 mais aucunement limitative a la determination d ! une 
quantite initiale d'acides nucleiques dans un echantillon 
soumis a une reaction de PCR (pour "Polymerase Chain 
Reaction"), en temps reel. Une technique de ce type, dite 
"de quantification par PCR", est utilisee notamment pour 
15 evaluer le nombre de copies d' agents pathogenes (par 
exemple du Virus de 1' Immunodef icience Humaine VIH) dans 
un pr61evement de fluides corporels d'un patient, 
typiquement dans le cadre d ? un suivi medical. 

20 On se refere a la figure 1 pour decrire brievement 
• 1 ? allure schematique d ! une courbe d 1 amplification par PCR 
en temps reel, avec, en abscisses, des indices de cycles 
de PCR et, en ordonnees, dans 1 ! exemple represents, des 
quantitfes de fluorescence emises (en unite arbitraire) et 

25 mesurees a chaque cycle de PCR. On indique qu'a chaque 
cycle de PCR, 1' echantillon subit des variations de 
temperature permettant a 1'ADN polymerase de proceder a 
1' amplification des acides nucleiques et les produits PCR 
correspondants sont detectes par des molecules 

30 f luorescentes . On obtient, en reportant la quantite de 
fluorescence mesuree F n en fonction du nombre n de cycle 



1er depot 



2 

de PCR, une variation du type reprdsente sur la figure 1 
et comprenant au moins : 

- une premiere partie BDF ou les mesures de fluorescence 
se confondent sensiblement avec un bruit de fond de 
l'appareil de mesure de fluorescence, 

- une seconde partie EXP ou les quantites de fluorescence 
mesurees croissenl de fagon sensiblement exponentielle, 

- une troisieme partie LIN ou la croissance des quantites 
de fluorescence mesurees est sensiblement attenuee et 
se comporte globalement de fagon sensiblement. lineaire, 
et 

- une quatrieme partie PLA ou les mesures de fluorescence . 
atteignent une phase de plateau. 

On retiendra en fait que pour les premiers cycles de PCR 
(premiere et seconde parties), la population d'interet 
croit sensiblement de fagon exponentielle, tandis que pour 
les cycles qui suivent (troisieme et quatrieme parties), 
d'autres phenomenes entrent en competition avec la 
croissance de la population d'interet, cette croissance 
etant alors amortie jusqu'a la phase de plateau PLA. 

On connait notamment par le document : 

"Mathematics of quantitative kinetic PCR and the 
application of standard curves", de R.G.Rutledge et 
C.Cote, dans Nucleic Acids Research , 2003, vol.31, N°16, 
une methode pour estimer la quantity initiale, inconnue, 
d f acides nucleiques dans un echantillon d'interet, par 
PCR. Elle consiste a utiliser plusieurs echantillons de 
quantites initiales connues en acides nucleiques, dits 
"standards", pour determiner, par interpolation, la 
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quantity initiale d'acides nucleiques presents dans 
1 'echantillon d'interet. 

De fa?on generale, plus la quantite initiale en acides 
5 nucleiques dans un echantillon est elevee, et plus tot la 
quantite de produit PCR obtenue et, par consequent, la 
fluorescence emise deyiendront d6tectables. En se r6f§rant 
a la figure 2 relative a I'art anterieur, on comprendra 
que la population initiale dans le standard Stl est plus 
10 elevee que celle du standard St2, laquelle est plus elevee 
encore que celle du standard St3, etc, car le cycle Ctl 
pour le standard Stl intervient avant le cycle Ct2 
correspondant pour le standard St2, lequel intervient 
encore avant le cycle Ct3 du standard St3, etc. 

15 

Ainsi, ce cycle Ct/< correspondant au cycle auquel la 
mesure de fluorescence atteint un seuil de fluorescence 
THR (comme repr^sente sur la figure 2) , prefixe 
arbitrairement ( typiquement au-dessus du bruit de fond), 

20 se presente comme un param£tre representatif de l'effectif 
initial N 0 d f une population d'acides nucleiques soumise 
aux cycles de PCR. Cette observation a ete mise a profit 
dans l'art anterieur cite ci-avant pour etablir une 
dependance telle que representee sur la figure 3 entre les 

25 numeros de cycle Ctl, Ct2, Ct3, Ct4, pour plusieurs 
standards d'effectifs initiaux connus, et leur effectif 
initial No 1 , N 0 2 , N 0 3 , N 0 4 . En reportant ainsi. les cycles 
Ctl, Ct2, Ct3, Ct4, etc, en ordonnees et le logarithme des 
effectifs initiaux N 0 X , N 0 2 , N 0 3 , Nq 4 en abscisses, on 

30 obtient une pente de regression REG. Sur cette pente REG, 
on reporte (fleche Fl en traits pointilles) le cycle Ctint 
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detecte pour 1' echantillon d'interet. Par interpolation 
sur la pente de regression REG (fleche F2 en traits 
pointings), on determine ensuite I'effecrif initial N 0 int 
de la population dans 1' echantillon d'interet. 

Cette methode, merne si elle s'est largement r<§pandue, 
pr6sente neanmoins quelques inconvenients . 

Tout d'abord, elle necessite 1 ' utilisation de plusieurs 
echantillons standard d'effectifs initiaux respectifs 
connus . 

Ensuite, la methode depend du juge de 1' utilisateur car la 
valeur du seuil de fluorescence qui est choisie par 
1' utilisateur influence directement la valeur des cycles 
Ct des courbes d' amplification et, par consequent, la 
valeur estimee de l'effectif initial dans 1' echantillon 
•d'interet. La valeur du seuil a egalement un impact sur la 
precision du resultat car la precision est generalement 
meilleure si le seuil est choisi dans la phase EXP cie 
croissance exponentielle des courbes d' amplification . 
Toutefois, en pratique, il est difficile pour 
1' utilisateur de savoir si le niveau seuil de fluorescence 
fixe THR correspond bien a la phase exponentielle des 
courbes, et ce, . pour tous les echantillons (les 
echantillons standard et 1 1 echantillon d f interet). 

Enfin, la methode sous-entend, sans aucune verification, 
que la population a le meme rendement d' amplification dans 
1' echantillon d'interet et dans tous les echantillons 
standard. Ainsi, si typiquement 1' echantillon d'interet 
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contient des inhibiteurs de PCR, son resultat sera 
f aussement abaisse . 

On 'retiendra done que la technique ' de l'art ant6rieur 
5 depend du seuil de fluorescence THR defini par 
l'utilisateur . La valeur choisie influence les valeurs des 
cycles Ct et par consequent la determination de la 
quantite initiale dans 1 » echantillon d'interet. C'est 
l'une des raisons pour lesquelles de nombreux travaux ont 
10 cherche recemment a automatiser et a fiabiliser la 
detection du cycle Ct . 

La prdsente invention vient ameliorer la situation en 
proposant une approche completement differente. 

15 

Elle vise tout d ' abord un procede, ce procede etant mis en 
ceuvre par des moyens inf ormatiques pour quantifier, de 
fagon absolue et/ou relative, une population initiale 
d'acides nucleiques dans un echantillon d'interet. Cet 

20 echantillon est soumis a une succession d 1 applications 
d J une reaction d 1 amplification de la population d ! interet. 
De fagon tres generale, cette amplification peut etre 
menee par cycles successifs cle PCR, mais aussi par toute 
autre technique d 1 amplification . On retiendra surtout que 

25 cette amplification doit simplement etre definie par un 
rendement de reaction, comme on le verra plus loin. Au 
cours de ces amplifications successives, on releve des 
mesures experimentales representatives d'un effectif 
courant de la population au moins dans 1 ' echantillon 

30 d ! interet. On comprendra que 1 ' on peut relever une ou 
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plusieurs mesures apres ou pendant chaque reaction 
d* amplification, sans perte de generality. 

Selon une definition actuellement prefer£e de 1' invention, 
le procede au sens de 1' invention comporte les etapes 
suivantes : 

a) prevoir un modele de rendement de la reaction 
d' amplification en fonction de la succession 
d T amplifications, ce modele comportant : 

o une phase sensiblement constante pour une premiere 

partie des applications de la reaction, et 
o une phase non constante pour une seconde partie 
des applications de la reaction, 
les premiere et seconde parties 6tant separees par une. 
region de basculement du rendement entre les phases 
constante et non constante, incluant un indice 
d' application correspondant sensiblement au 

basculement, 

b) exprimer, en utilisant le modele de rendement, une. 
relation faisant intervenir au moins 1' indice de 
basculement et un parametre representant 1 ' ef f ectif 
initial dans 1 ' echantillon d'inter§t, 

c) par comparaison avec les mesures experimentales, 
determiner au moins 1? indice de basculement, et, dans 
une etape d) ulterieure ou conjointe, en d£duire 
l T effectif initial dans 1 1 echantillon d'interet. 

D ! autres avantages et caracteristiques de 1' invention 
apparaitront a la lecture de . la description detaillee du 
mode de realisation donne ci-apres a titre d f exemple, en 
s'appuyant sur les graphiques annexes et sur lesquels : 
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outre la figure 1, relative a l'art anterieur, 
representant une variation de la quantity de 
fluorescence- mesuree en fonction d'un nombre de cycles 
de PGR, decrite ci-avant ; 

outre la figure 2, relative a l'art anterieur, 
representant les quantites de fluorescence <§mises et 
croissantes en fonction du nombre de cycles de PCR, et 
decrite ci-avant ; 

outre la figure. 3 representant une m<§thode 
d 1 interpolation pour determiner la quantite initiale de 
la population d'interet dans 1 1 echantillon d'interet, 
selon une methode connue de l'art anterieur, et decrite 
ci-avant ; 

la figure 4A represente une variation schematique des 
mesures exp^rimentales precitees en fonction de la 
succession d 1 amplifications que 1 1 on applique & 
1 1 echantillon d'interet, 

la figure 4B represente une. variation schematique du 
rendement de la reaction d ' amplification, obtenue a 
partir des mesures experimentales , en fonction de la 
succession d ' amplifications appliquees a 1 1 echantillon 
d' interet , 

la figure 5 represente une regression entre des indices 
de basculement du rendement entre la phase constante et 
la phase non constante et le logarithme des effectifs 
initiaux pour des echantillons standard et pour 
1' echantillon d'interet, pour la. mise en oeuvre d r un 
premier mode de realisation, 

la figure 6A represente une variation typique de la 
quantite de fluorescence mesuree et ajustee en tenant 
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compte d'un bruit de fond propre aux mesures, en 
fonction du nombre de cycles n de la reaction PCR, 

- la figure 6B represente une variation de 1 1 ef f icacite 
de la PCR de la figure 6A, en fonction du nombre de 
cycles n, 

- la figure 7 represente une comparaison de la variation 
experimentale de la fluorescence emise selon la figure 
6A et le resultat obtenu par 1 ' application du modele de 
la fluorescence emise obtenu par integration d ! un 
modele d 1 eff icacite, selon un second mode de 
realisation, 

- la figure 8 represente une comparaison de la variation 
de l'efficacite de la figure 6B avec 1 1 application d'un 
modele d 1 eff icacite, deduit de la modelisation de la 
fluorescence emise de la figure 7, 

- la figure 9 represente 1 ' organigramme schematisant les 
principales etapes du procede selon le mode de 
realisation propose pour la presente invention, et 

- la figure 10 represente une installation pour 
quantifier la population initiale d'un echantillon 
d'interet. 

On s'appuie sur les figures 4A et 4B pour decrire 
brievement quelques principes de 1' invention illustrant 
les caracteristiques du procede ci-avant. 

Tout d'abord on indique que la figure 4A represente une 
succession de mesures experimentales F n representatives de 
l'effectif courant d f une population d f int6ret qui subit 
progressivement une succession de reactions 

d ! amplif ication, chaque reaction etant indexee par un 
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indice n. Cette succession de reactions correspond, dans 
l f exemple decrit ici a titre non limitatif, a une 
succession de cycles de PCR. Les mesures experimentales F n 
correspondent dans cet exemple a titre non limitatif a des 

5 mesures de quantites de fluorescence a chaque cycle de 
PCR. Ainsi,. dans le cadre d'une methode de quantification 
qui* combine la reaction PGR et la fluorescence ernise par 
1 1 echantillon d'interet, on introduit des reactifs 
fluorescents dans 1 ' echantillon, de sorte que la 

10 fluorescence emise, au cours d'un cycle de PCR, est 
proportionnelle a I'effectif d'acides nucleiques dans cet 
echantillon. . Neanmoins, cette realisation admet de 
nombreuses variantes. Deja, il peut etre pref6re de mener 
plusieurs mesures ou aucune mesure pour certains cycles de 

15 • PCR. Par ailleurs, plus generalement , la methode de mesure 
peut utiliser d'autres techniques que la fluorescence, 
meme si cette derniere est souv.ent celle employee dans la 
quantification par PCR. Enfin, on indique que d'autres 
techniques d 1 amplification peuvent etre prevues dans le 

20 cadre de la presente invention, d&s lors qu'il est 
possible de suivre la variation du rendement de la 
reaction correspondant a cette amplification. Neanmoins, 
comme l 1 exemple decrit ci-apres traite pref erentiellement 
de cycles de PCR, on parlera d ! " efficacite de la PCR", 

25 not^e E n pour chaque cycle de PCR d' indice n, pour 
designer le rendement de la reaction d 1 amplification . 

Comme indique ci-avant en reference a la figure 1, on 
retrouve sur la figure 4A principalement deux regions ou : 
30 - pendant les premiers cycles de PCR (partie EXP) , la 
population croit sensiblement de fagon exponentielle, 
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- tandis que pour les cycles qui suivent (partie LIN et 
PLA), d ! autres phenomenes entrent en competition avec 
la croissance de la population d f interest, cette 
croissance etant alors amortie. 

On formule alors les hypotheses suivantes : 

- le rendement de la reaction E n est relat ivement constant 
durant les premiers cycles sur la partie EXP, et . 

- a partir d f un certain nombre de cycles effectues, le 
rendement E n de la reaction decroit sur la partie LIN et 
PLA . 

Cette decroissance du rendement peut avoir diverses 
explications, notamment une degradation et/ou une carence 
des r6actifs de PGR (DNA polymerase, dNTPs, amorces, ou 
autres) et/ou une inhibition par les produits formes eux- 
memes. 

On fait ainsi l'hypothese ici que le rendement est d'abord 
constant et decroit ensuite . . Toutef ois, on comprendra que 
1' invention s 1 applique de facon plus generale dans de. 
cadre d'un rendement : 

- d'abord constant, ce qui correspond k une situation 
normale de croissance par amplification, 

- puis non constant (decroissant ou croissant), ce qui 
correspond a une situation sensiblement anormale. 

On a constate que dans le cadre de reactions 
d' amplification de quantite d'acides nucleiques, le 
rendement bascule souvent d'une phase constante a une 
phase non constante. Cette observation a ete mise a profit 
au sens de 1 ! invention pour en ■ deduire la quantite 
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initiale d'acides nucleiques, comme on le verra en detail 
ci-apres. Auparavant, ■ on indique simplement que le 
rendement peut encore basculer d ! une phase non constante 
pour les premiers cycles a une phase constante ensuite. La' 
5 presente invention s' applique tout aussi bien a ce cas . De 
maniere generale, on retiendra done que 1 ' on detecte, au 
sens de l 1 invention, un basculement du rendemenL, d'une 
phase constante a une phase non constante. 

10 On rappelle que le but recherche est de donner 1' effectif 
initial de la population soumise a amplification. Deja, en 
se ref£rant a la figure 4A, on comprendra que la mesure F 0 
representative de cet effectif initial, pratiquement 
confondue avec le bruit de fond des mesures BDF ne permet 

15 pas directement d'acceder, a elle seule, a l f effectif 
initial. Dans l'art anterieur, on a recherche a quantifier 
cet effectif initial en exploitant la- phase exponentielle, 
done typiquement a la sortie du bruit de fond. On y 
determine alors le cycle seuil Ct (correspondant au point 

20 AA pour "art anterieur n ) de la figure 4A. Or, comme on l'a 
vu ci-avant, les mesures sont souvent bruitees dans cette 
region et il est difficile de determiner avec precision un 
cycle seuil Ct de sortie du bruit de fond. 

25 Dans une toute autre approche, la presente invention 
exploite plutot la quasi-totalite des points de la courbe 
d ' amplification pour determiner precis6ment une region BAS 
oil le rendement bascule entre une phase constante et une 
phase non constante, typiquement ici entre la phase 

30 exponentielle EXP et la phase lin6aire LIN. Deja, on 
constatera que les mesures sont logiquement moins bruitees 
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dans cette region BAS que dans la region de sortie du 
bruit de fond puisque la region BAS intervient pour des 
cycles, ulterieurs. De plus, notamment en raison des 
proprietes mathematiques liees au rendement, on verra plus 
loin que le nombre de standards qui peuveht etre utilises 
pour la quantification de l'effectif initial dans la 
population d'interet est, de fagon tres avantageuse, 
reduit par rapport' au nombre de standards utilises dans la 
quantification de l'art anterieur. 

On decrit maintenant brievement la relation qui peut lier 
la region de basculement BAS a l ! effectif initial de la 
population d'interet. On rappelle au prealable que le 
rendement d ! une reaction d 1 amplification est donne par : 

N n+ 1 - N n + E n .N n , OU I 

- N n est l'effectif de la population d'interet apres une 
amplification d'indice n dans une succession 
d' amplifications, 

- N n+ i est l'effectif de la population d'interet apres une 
amplification suivante, d'indice n+1, dans la 
succession d ' amplifications precitee, et 

- E n est le rendement de la reaction d 1 amplification 
d'indice n dans la succession d 1 amplifications 
precitee . 

En reformulant cette relation, par recurrence, on 
obtient : 

Nn+1 = (1 + E n ) (1 + E n -i) (1 + E„- 2 ) ... (1 + E 0 ) N 0 , 
ou N 0 est l'effectif initial de la population d'interet. 
Tant que le rendement E n est constant, on comprendra que 
cette derniere relation peut s'ecrire simplement : 
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N n+ i - (1 + E 0 ) n+1 . N 0 (ou l'indice n+1 n'a pas 
encore atteint la region de basculement BAS) . 
En effet, le rendement etant alors constant pour ces 
premiers cycles, on a : 

En ~ E n -i = E n -2 — ... ~ E(), 

ou E 0 est la valeur de la phase constante du rendement. 
Neanrnoins, lorsque l'indice n+1 d6passe la' region de 
basculement BAS, la relation devient : 

N n +i = N 0 x (1 + E 0 ) C eep x fonction (C FfEP , n+1), ou : 

- (C EE p - 1) est le dernier indice de la reaction 
d 1 amplification au cours duquel le rendement est encore 
constant (on comprendra done que l'indice C EE p 
represente lui-meme l'indice de basculement, proprement 
dit, entre la phase de exponentielle a la phase 
lineaire ) , 

- le terme fonction (C EEP , n+1) est une fonction 
particuliere qui caracterise la phase non constante du 
rendement et qui depend au moins de l'indice de 
basculement C EEP et de l'indice courant de 
1 1 amplification n+1. 

On, voit ainsi comment peuvent etre lies l'indice de 
basculement C EE p et . 1 1 ef f ectif initial N 0 de la population 
d'interet. On comprendra qu ' a ce stade, on a deja realise 
les etapes a) et b) du procede tel que defini ci-avant. 

Deja, un premier mode de realisation consiste a determiner 
experimentalement l'indice de basculement C EEP et a • le 
correler a l'eff ectif initial, par regression, en 
utilisant plusieurs echantillons standard qui subissent le 
meme traitement - d ' amplification que 1 1 echantillon 
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d'interet. Dans ce cas, on comprendra que les etapes b) et 
c) du procede tel que defini ci-avant sont simplement 
inversees puisque l'on determine experimentalement , 
d'abord, l'indice de basculement C EEP (etape c)) et que 
l f on determine ensuite la relation entre l'indice C EEP et 
I'effectif initial N 0 (etape b) ) pour obtenir finalement 
I'effectif initial N 0 (6tape'd)). 

Avant de decrire en detail toutes ces etapes au sens de ce 
premier mode de realisation, on decrit un procede .de 
determination de l'indice C EEP , a partir des mesures 
experimentales. En particulier, on comprendra que cette 
maniere de determiner experimentalement l'indice C EE p peut 
etre appliquee a un autre mode de realisation, different 
du premier mode de realisation precite. 

En reprenant la relation entre l'efficacire E n a un cycle 
donne n et I'effectif courant dans 1 1 6chantillon d'interet 
au meme cycle N n et a un cycle ulterieur N n+ i, on peut 
exprimer l ! efficacite de 1 1 amplification comme suit : 
E n = N n+ i/N n - 1 

Dans certains cas. notamment ou il n'est pas necessaire de 
tenir compte du bruit de fond BDF dans les mesures, on 
peut en premiere approximation considerer que les mesures 
sont sensiblement proportionnelles a I'effectif courant 
dans l'eqhantillon d'interet. Toutefois, en pratique, on 
tiendra le plus souvent compte des derives des mesures et, 
dans une realisation qui sera decrite maintenant, on 
determine des mesures experimentales corrigees, notees 
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F ! n , a partir ■ des mesures directes F n telles que 
representees sur la figure 4A. 

Pref erentiellement, on applique une etape prealable de 
5 traitement des mesures experimentales F n consistant a 
soustraire le bruit de fond BDF et £ introduire ensuite 
une compensation pour tenir compte d'une mesure non nulle 
e representative de l ! effectif initial, Dans l'exemple 
represents sur la figure 4A, la variation du bruit de fond 

10 BDF en fonction de 1-' indice n peut etre representee par 
une fonction lineaire car les essais ont montre qu'un 
modele lineaire est satisfaisant pour des mesures de 
fluorescence en PCR. N<§anmoins, un modele de variation 
exponentielle pourra etre pref6re dans certains cas. On 

15 applique en tout ■ etat de cause le modele qui respecte le 
mieux la variation du bruit de fond BDF donnee typiquement 
par les premiers points de mesure. Ensuite, on soustrait 
le modele BDF de variation de bruit de fond choisi a 
toutes les valeurs de mesures experimentales ■ F„. En 

20 appliquant cette etape, on comprendra que la mesure de 
fluorescence theorique Fq est ramenee a une valeur de 
mesure nulle, ce qui correspond a un effectif initial N 0 
nul, ce qui n'a aucune realite physique. Par consequent, 
on applique avantageusement une compensation de . cette 

25 correction comme suit : 

F ! n = F n - BDF + e, ou : 

- le terme F f n correspond a une mesure corrigee pour un 
indice courant n, 

- le terme F n correspond a la mesure experimentale brute a 
30 cet indice courant n, 
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- le terme BDF correspond a la valeur du bruit de fond 
modelise pour l'indice courant n, et 

- e est le terme de compensation correspondant , suppose 
constant, dans l'exemple decrit et representant 
directement en fait l'effectif initial N 0 . 

Ces etapes de correction par rapport au bruit de fond, 
meme si elles sont tres avantageuses pour la determination 
de l'indice de basculement C EE p, peuvent etre appliquees 
par ailleurs a toute determination et quantification de 
l'effectif initial N 0 des lors qu'un bruit de fond est 
susceptible de fausser la mesure de cet effectif N 0 . A ce . 
titre, ces etapes peuvent faire l'objet d'une protection 
separee, le cas echeant. 

Les mesures corrigees F' n ainsi obtenues. sont 
avantageusement proportionnelles auk effectifs cour.ants N n - 
dans 1 ' 6chantillon d'int^ret, de sorte que le rendement E n 
peut maintenant s'exprimer directement en fonction de.s 
valeurs de mesure (corrigees comme indiqu£ ci-avant) , par 
la relation : 

E n - F'n^/F'n - 1 

Ainsi, a partir des mesures experimentales F n de la figure 
4A, on obtient les mesures corrigees F' n a partir 
desquelles ' on determine ensuite la variation de 
1'efficacite E n en fonction de la succession d f indices n, 
telle que representee sur la figure 4B. 

En bref, on exprime les mesures experimentales sous la 
forme d ! une variation experimentale du rendement E n telle 
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que representee sur la figure 4B, en fonction de la 
succession d 1 applications n. On obtient alors une 
variation experimental, du rendement, comportant : 

- une premiere region sensiblement bruitee pour des 
5 indices d 1 applications n faibles (sensiblement avant le 

cycle CG dans l'exemple de la figure 4B) , 
--suivi.e d T une seconde region moins -bruitee pour des 
indices d ' applications plus eleves (au moins apres la 
region de basculement BAS) .■ 

10 

Au moins dans le cas le plus courant en amplification par 
PCR et en mesure par fluorescence, la phase non constante 
du rendement est decroissante et correspond a cette 
seconde region peu bruitee (comme represents sur la figure 
15 4B) . Notamment pour eliminer des points de mesure qui 
risqueraient de fausser lea resultats lorsque 1 1 on choisit 
un module a appliquer a la variation du rendement : 

- on estime une valeur grossiere de la phase constante du 
rendement E 0 , et 

20 - notamment pour la recherche de 1 ' indice de basculement 
C E ep, on ignore une partie au moins des mesures dans la 
seconde region moins bruitee et dont le rendement 
estime est inferieur £ une valeur seuil, telle que, par 
exemple, une fraction de la phase constante E 0 . 

25 

Ces points NEG (figure- 4B) qui sont elimines sont 
typiquement ceux correspondant a des indices 
d 1 amplification n tres eleves et qui pourraient ne plus 
respecter le modele qui a ete choisi pour l'efficacite 
30 sensiblement autour de la region de basculement BAS. Pour 
les eliminer, on evaluera" par exemple une moyenne de la 
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phase constante du rendement E 0 typiquement pour les 
premiers indices n. Ensuite, on choisira une valeur seuil 
correspondant a une fraction, par exemple 10%, de la 
moyenne trouvee pour la phase constante E 0 . Ensuite, en 
partant des indices n les plus eleves, on §limine tous les 
points de mesure NEG dont la mesure de rendement est 
inferieure ou egale a cette valeur seuil. Cette etape, 
tres avantageuse pour la detection de 1 1 indice Geep* peut 
nfeanmoins etre appliquee a toute determination basee sur 
le rendement E n et, le cas echeant, peut faire l'obj.et 
d'une protection separee. 

Dans le" cas d'une phase non constante du rendement qui est. 
decroissante et qui suit la phase constante, comme 
represents sur la figure 4B, on repere la region de.. 
basculement BAS, dans le sens des indices n decroissant s 3 
partir de la seconde region moins bruitee, en detectant un r 
indice grossier CG dont le rendement depasse une valeur 
predeterminee. Ainsi, en se referant a la figure 4B, on va 
dans le sens des indices n decroissants en remontant vers 
la region de basculement BAS et on evalue le rendement 
associe a cet indice. Pour le premier point de mesure dont 
le rendement depasse signif icativement la valeur 
predeterminee precitee, on considere que 1' indice grossier 
precite CG est detecte et correspond a 1 ! indice de ce 
point de mesure. 

Comme on le verra dans un mode de realisation ulterieur, 
on peut modeliser, pour chaque point de mesure, la 
variation de son rendement comme si ce point correspondait. 
effectivement a 1' indice de basculement C EE p- Dans cette 
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realisation, on estime alors sa phase constante E 0 de 
rendement et si cette valeur estimee depasse la . valeur 
predeterminee precitee, on considere ce point comme 
correspondant a l'indice grossier CG. 

5 

De fagon generale, un rendement maximum a une valeur de 1 
de sorte que i 1 on peut 1 choisir la valeur predeterminee 
precitee comme etant egale a 1. Neanmoins, ce choix admet 
des variantes et par exemple, il peut etre prevu de fixer 
io cette valeur predeterminee comme correspondant, a la 
moyenne du rendement E 0 evaluee sur les premiers cycles de 
reaction . 

Ensuite, on affine l f estimation de la valeur, qui peut 

15 avantageusement etre f ractionnaire, ' de l'indice 
d T application C E ep dans la region de basculement, dans le 
sens des indices d 1 applications croissants a partir de 
l'indice grossier CG, en detectant un indice d 1 application 
d' amplification dont le rendement vaut environ la valeur 

20 predeterminee precitee . Ainsi, en se referant a nouveau a 
'la figure 4B, pour af finer la recherche de l'indice de 
basculement C EEP apres la determination de l'indice 
grossier CG, on redescend a partir de l'indice grossier CG 
dans le sens des indices n croissants, par pas plus petit 

25 qu'un indice entier et on determine, par exemple par 
interpolation, l'abscisse pour laquelle on franchit la 
valeur predeterminee. Typiquement, tant cue la valeur de 
phase constante reste superieure a 1, on poursuit la 
recherche dans le sens des n croissants et 1'on determine 

30 l'indice (C EEP - 1) qui precede le basculement des que la 
valeur constante E 0 est egale ou est tres proche de 1. 
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C'est la raison pour laquelle on. choisira un pas de 
recherche correspondant a une fraction, d' indice, par 
exemple 10% d'un cycle n. 

Dans le premier mode de realisation precit6, on prevoit 
une pluralite d 1 echantillons standard d'effectifs initiaux 
respectifs connus, auxquels on applique la" succession 
d' amplifications, sensiblement dans les memes conditions 
que pour 1 ' echantillon d'int6ret. On determine aussi leurs 
indices respectifs de basculement selon les Stapes a) , b) 
et c) decrites precedemment , A 1 ' etape d) : 

- on etablit une dependance entre les effectifs initiaux 
des echantillons standard N 0 st et leurs indices C EEP st , 

- et, apres la determination de 1 ' indice pour 
• l» echantillon d'interet C EEP , on determine l'effectif 

initial de la population d'interet N 0 par interpolation 
sur cette dependance. 

Ainsi, en se r^ferant a la figure 5, on peut etablir une 
dependance entre le cycle de basculement Ceep^Ceep 2 des 
standards et leur concentration initiale No^No 2 (en fait 
leur logarithme), par exemple par regression. En mesurant 
1» indice de basculement C EE p int pour 1 1 echantillon d ! interet 
et en reportant sa valeur sur la pente de regression de la 
figure 5, on obtient par interpolation la concentration 
initiale N 0 int dans 1 ' echantillon d'interet. 

Ce premier mode de realisation s'apparen.te done 
sensiblement a celui de l'art anterieur decrit en 
reference a la figure .3. Cependant, on gardera a 1' esprit 
que 1 ! indice de basculement C EEP sur lequel se base ce 
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premier' mode de • realisation ne correspond aucunement au 
cycle Ct de l f art anterieur. 

Dans une approche sensiblement differente de ce premier 
mode de realisation : 

- a l'etape b) , on exprime en fait, en utilisant le 
■ modele de 'rendement, une variation parametree en 

fonction de la succession d' amplifications, cette 
variation mettant en jeu au moins un parametre 
representant l'indice de basculement C EEP , et 

- a l'6tape c) , par comparaison avec les mesures 
experimentales, on determine au moins ce ' parametre 
representant l'indice de basculement C EEP . 

Dans un second mode de realisation, cette variation 
parametree est representative de l ! effectif courant N n 
dans 1 r echantillon d*interet. 

Typiquement, cette variation parametree peut etre tiree de 
1' expression du type donne ci-avant : 

N n+1 = N 0 x (1 + E 0 ) C eep x fonction (C EEP , n+1) 
Ainsi, cette variation met en jeu, outre un parametre 
representant l'indice de basculement C EE p, un parametre 
representatif de l'effectif initial dans 1 1 echantillon 
d'interet N 0 - 

Ensuite, aux etapes c) et d) de ce second mode de 
realisation, on determine sensiblement conjointement ces 
deux parametres C EEP et N 0 . 
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Auparavant, a l'etape a), il faut determiner un modele 
pour la fonction precitee fonction (C E ep, n+1) . 

Le plus souvent, pour une quantification par PCR, on 
choisira un modele pour la phase non constante du 
rendement correspondant a une exponentielle decroissante 
comportant un parametre de ■ decroissance 3 qui sera plus 
detaille ci-apres. Ce parametre de decroissance p est 
alors determine, a l'etape c) , au moins avec 1 ' indice de 
basculement C EE p, P*r comparaison avec les mesures 
experimentales . 

Ainsi, dans ce second mode de realisation, une fois le 
modele du rendement E n choisi, on 1' applique a 
1' expression generale de l'effectif courant N n donnee par 
la relation ci-avant. On obtient alors un modele de la 
variation de l'effectif courant N n . 

Toutefois, a moins que les mesures experimentales ne 
donnent directement la valeur de l'effectif courant N n (ce 
qui est rarement le cas en pratique, actuellement ) , il 
convient de modeliser ensuite les mesures experimentales 
F n , elles-memes, en gardant a l 7 esprit la soustraction du 
bruit de fond et la compensation consecutive e decrites 
ci-avant . 

Ainsi, dans une realisation actuellement preferee/ la 
variation parametree precitee : 

- est representative des mesures experimentales, et 

- comporte un parametre correspondant- a une valeur de 
mesure F 0 representative de l'effectif initial. 
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On determine ensuite la valeur de mesure de l ! effectif 
initial F 0 par comparaison de cette variation parametree 
F n avec les mesures experimentales. 

En fait, pour mener cette comparaison, on pourra par 
exemple ajuster les parametres F 0/ E 0/ C EEP et le parametre 
de decroissance p dans le modele des mesures F n par 
correlations statistiques (typiquement la m^thode des 
moindres carres) qui s'appuient sur les mesures 
experimentales brutes. Un exemple de realisation est 
d§ta"ille ci-apres. 

On obtient d'abord une variation de la quantite de 
fluorescence mesuree et ajustee, en fonction du nombre de 
cycles de PCR appliques, telle que representee par exemple 
sur la figure 6A. Cette figure montre la courbe 
d f amplification d'un echantillon d'int6ret contenant des 
acides nucleiques, ici un fragment d'ADN de quantite 
initiale de 100000 copies, marque par 1 ' intercalant 
SYBRGREEN lors de la reaction de PCR que 1 ■ on effectue sur 
l'appareil I-Cycler IQ® de la society BIO-RAD®. 

Dans 1' exemple decrit, on comprendra que la reaction 
d ! amplification est une reaction de PCR en temps reel. Les 
mesures experimentales representent des quantites de 
fluorescence emises . 

On note ci-apres F n la fluorescence au cycle n 
prealablement ajustee pour le bruit de fond, comme indique 
ci-avant. On note aussi F 0 la fluorescence initiale 
theorique avant le premier cycle. On note E n I'efficacite 
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de la PCR au cycle n. On d6signe par N le nombre total de 
cycles effectues durant la reaction de PCR*. 

Par hypothese, la fluorescence mesuree a chaque cycle n de 
la reaction de PCR est definie par : 

F n+l * F n (l+E n ) pour tout ne{0,l,2,...,N-l} (1) 
avec 0 ^ E n ^ 1* 

On calcule l'efficacite de la reaction a chaque cycle n, 
comme suit : 

E = - 1 « ^ - 1 pour tout n e {l,2,...,N-l} ( 2 ) 

N F 

On remarquera que 1' equation (1) est supposee vraie pour 
n=0. Toutefois, par definition, la fluorescence initiale 
F 0 est inconnue. L'efficacite au premier cycle E 0 ne peut 
done pas §tre calculee directement par la formule (2) . 

On a represents sur la figure 6B l'efficacite de la 
reaction PCR approximee par la formule (2) et a partir de 
la variation de la fluorescence ajustee de la figure 6A, 
en fonction du nombre de cycles n. 

On. formule pref erentiellement les hypotheses suivantes : 

- l'efficacite de la reaction est relativement constante 
durant les premiers cycles, et 

- a partir d'un certain nombre de cycles effectues, 
l'efficacite de la reaction decroit. 

La figure 6B vient conforter la seconde hypothese, puisque 
l'on peut constater la decroissance de l'efficacite a 
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partir du cycle n=17 . En revanche, les efficacites des 
cycles 1 a 16 sont fortement bruitees ce qui rend 
difficile la verification graphique de la premiere 
hypothese. 

Toutefois, on affecte pref erentiellement a la variation de 
l'efficacite un modele du type incluant : 

- une premiere phase, constante, entre le premier cycle 
de PCR et le cycle (C EEP -1) qui precede le cycle de 
basculement note C E ep/ et 

- une seconde phase, decroissante , pour les cycles 
superieurs ou egaux au cycle (C EEP -1) . 

Le cycle (C EEP -1) represente done le dernier cycle, qui 
peut etre f ractionnaire, auquel l'efficacite est encore 
constante . 

On propose alors de modeliser l'efficacite de la reaction 
de la fagon suivante : 



ou E 0 et p sont des parametres reels qui sont estimes 
grace a la courbe d 1 amplification de la figure 6A, ou 
encore grace a la courbe d'efficacite de la figure 6B, 
comme on le verra ci-apres. 

En variante, on indique qu ! un autre choix, par exemple 
parmi les modeles ci-apres Fl a F3, peut etre prefere, 
notamment en fonction du type d'acide nucleique que l f on 
souhaite quantifier : 




pour O^n <(C EEP -1) 
pour (C EEP ~l)<n <; (N-l) 



(3) 
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Fl : 


E„ 


= exp(-p(n-C EEP +l))-l 


F2 : 


E„ 


= exp(-|i(n-C EEP +l) a )- 1 


F3 : 


E n 


= a-exp(-n(n-C EEP +l) B )- 1 



On estime pref erentiellement plusieurs jeux de parametres 
a 1 ' etape c) pour plusieurs candidats de cycle de 
5 basculement C EEP , et 1 } on choisit le cycle candidat minimal 
dont les parametres associes maximisent les correlations 
statistiques qui peuvent etre menees a l f £tape c) , en tant 
que cycle de basculement C EEP . 

10 Comme on l'a vu precedemment , 1' expression (1) • ci-avant 
peut encore s'ecrire sous la forme : 

F„ = Fo flO+Ej Pourn e {l,2,..,N} (4) 

Ainsi, en introduisant 1 ! expression de l'efficacite (3) a. 
15 dans la formule (4), on obtient un' nouveau modele a quatreg 
parametres (F 0 , *E 0 , P, C EEP ) pour la fluorescence emise" 
ajustee l n • 

f F 0 (l+E 0 ) n pour l<n < C EEP 

F 0 (l+E 0 ) UKV - P (p)-i pour C EEP < n < N 

20 L'effectif initial dans la. population d'interet N 0 , 
l'efficacite E 0 de la reaction a n=0, le parametre (3 et le 
cycle de basculement C BE p, sont evalues repetitivement pour 
plusieurs valeurs de cycles dans la region de basculement 
BAS, pour trouver un maximum de correlations statistiques 

25 atteint pour une valeur de cycle minimale egale au cycle 
de basculement C EEP . 
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Dans ce second mode de realisation, on prefere modeliser 
la variation des quantites de fluorescence mesurdes et 
ajustees, en fonction des cycles, en se basant sur la 
5 modelisation de la variation de 1 1 eff icacite, et mene'r 
ensuite les correlations directement sur les quantites de 
fluorescence mesurees et ajustees. 

On remarquera qu ? en ajustant les mesures de fluorescence 
10 eraises pour le bruit de fond, on cree aussi, de fagon 

factice, un ajustement de la fluorescence initiale F 0 . 

Ainsi, estimer les parametres du modele d" eff icacite £ 

partir de mesures d 1 eff icacite . qui sont deduites des 

mesures de fluorescence ajustees, est une source d'erreurs 
15 supplementaires et il peut etre prefere de proceder en 

deux etapes, cornme on le verra dans le troisieme mode de 

realisation qui sera decrit plus loin. 

Toutefois, le second mode de realisation propose est plus 
20 simple et s'adapte bien aux quantifications par PCR 
utilisant des mesures de fluorescence. II se base sur les 
mesures reelles de fluorescence F' n , qui correspondent aux 
mesures de fluorescence ajustees pour la derive du bruit 
de fond avec une compensation z men£e sur ces mesures. Le 
25 bruit de fond ayant ete soustrait, on retrouve une 
relation du type : 

ou e est une quantite qui peut dependre ou non du cycle n. 
Elle est choisie ici constant©, pref erentiellement . 
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Dans ce cas, l'eff icacite mesuree et "ajustee" aussi, 
notee E' n au cycle n est dfefinie par- : 

E 1 = ^ - 1 = Fn+1 + £ - 1 pour tout n e {l,2,...,N-l} ( 7 ) 

F' n F n + e 

Le modele de la relation (5) ci-avant devient : 
F' 0 (l+E' 0 ) n -e pour 1< 

F n = \ Lexp (-P '(n -C Ht!P )) ( 8 ^ 

F' 0 (l+E' 0 ) "* ctpivi - e pour C EEP < n <; N 

Dans ce cas., les valeurs d f efficacit6 E' n sont approximees 
experimentalement a partir des mesures, de fag:on a pouvoir 
fixer un seuil d'efficacite minimal toler£ dans la phase 
cTeff icacite dScroissante . On determine ainsi un cycle 
seuil au-dela duquel les mesures de fluorescence ajustees 
(points NEG figure 4B) ne seront pas exploiters pour mener 
la modelisation. Typiquement, le cycle seuil correspond au 
premier cycle de la phase d'eff icacite decroissante' pour 
lequel l'efficacite est inferieure au seuil d'efficacite 
minimal tolere (par exemple 0,1.E 0 ). 

Plus generalement, la valeur du seuil d'eff icacite peut 
pref erentiellement etre comprise entre 0 et 0,5, et une 
efficacite de PCR inferieure a cette valeur est 
potentiellement biaisee par des phenomenes d ! inhibition 
incontr616s. 

Dans 1' exemple represents sur la figure 7, la valeur seuil 
d'efficacite a ete fixee a 0,02 (soit 2% de E 0 ) . Le cycle 
seuil C s correspondrait au cycle n=36. On a reprSsente sur 
cette figure 7 les mesures de fluorescence emises 
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ajust6es. On constate une correlation satisf aisante du 
modele (trait plein) avec les mesures experimentales 
(marquees "o") jusqu'au cycle n=36. La figure 8 montre 
aussi une bonne correlation avec les mesures 
experimentales des predictions de 1 ' ef f icacite, obtenue a 
partir de la figure 7 et de la modelisation de la 
fluorescence mesurere et ajustee. 

Les etapes principales de ce mode de realisation peuvent 
se r£sumer comme suit, en reference a la figure 9. 



A l'etape de debut 70, on a obtenu les valeurs de 
quantites de fluorescence mesurees et ajustees par rapport 
au bruit de fond en fonction du nombre de cycles n, comme 
15 represents sur la figure 6A. 

A l'etape 71, on calcule une approximation de 1 1 eff icacite 
de la reaction au cycle n grace a la formule (2) ci-avant, 
pour tout cycle n = 1, 2, (N-l) . 



A I'etape 72, on determine le cycle C s minimal de fagon a 
ce que les deux conditions suivantes sont satisfaites : 

- le cycle C s se situe dans la phase d'efficacite 
decroissante, 

- 1' eff icacite au cycle seuil est inferieure a la valeur 
seuil d'efficacite E s (par exemple E s = 0,1. E 0 ) : 

On peut deja eliminer les points NEG pour lesquels 
l'efficacite est inferieure a E 3 . 
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A l'etape 73, on modelise la courbe de fluorescence 6mise 
ajustee dans la decroissance de 1 1 ef f icacite, sur la 
portion de cycles. C EEP = (C s - 5) a C s , grace a l'expression 
(8) dans laquelle on a fait l'hypothese sur la 
compensation e telle que s=F' 0 : 

l-exp(-PXn-C B Bp)) 

F n -F' 0 (l+E' 0 ) CREP -K»-i -F' 0 

Ensuite, le test 7 4 sur la valeur estimee E' 0 de E' 0 et la 
decrementation a l'etape 75 de la valeur du cycle de 
basculement C E ep visent a diminuer la valeur recherchee de 
C EEP d ! un pas P (qui peut etre egal a 1) et a repeter 
l'etape 73 tant que la valeur de E' 0 est inferieure a 1. 

Ensuite, lorsque la valeur d'efficacite estimee depasse la 
valeur 1 (fleche n en sortie du test 74), on incremente la 
valeur de l'indice C EE p d'un pas h (qui peut etre une: 
fraction d'unite) a l'etape 76 et on precede £ une etape; 
77 de modelisation de la fluorescence F n du type de celle 
menee & l'etape 73. Tant que l'eff icacite estimee E* 0 est 
super ieure ou egale a 1 au test 78 , on repete les e tapes 
76 a 78. Lorsque l'efficacite estimee prend une valeur 
inferieure a 1, on conserve les parametres estimes 

(FVE'oi 1 ^^) a l'etape de fin 79. 

A cette etape, on a obtenu finalement une valeur F' 0 
seulement representative de 1'effectif initial N 0 dans 
1' echantillon d'interet. On peut alors avoir recours a au 
moins un echantillon standard d'effectif initial N 0 st 
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connu, pour determiner, a I'etape 80, 1' effectif initial 
N 0 dans 1' echantillon d'interet. 

A cet effet, on obtient une valeur de mesure d ' un effectif 
initial F 0st dans un echantillon standard d f effectif 
initial connu N 0st . On determine ensuite la valeur de 
1' effectif initial dans 1 1 echantillon d ! int6ret N 0 par 
application d'une meme relation de proportionnalite entre 
l'effectif initial N 0 et sa mesure F' 0 , dans 1 ' Echantillon 
standard d ! une part, et dans 1 f echantillon d'interet 
d* autre part. 

En d'autres termes, a l'etape 80 de la figure 9, on peut 
determiner la valeur N 0 de 1' effectif initial dans 
1 1 echantillon d'interet par une simple regie de trois du 
type : N 0 = F 0 (N 0st /F 0sl ) , impliquant 1' effectif initial dans 

le standard N 0s t et un rapport des fluorescences corrigees, 
compensees et estimees par ajustement du module de 
fluorescence, pour 1' echantillon d'interet et pour 
1' echantillon standard. 

On • comprendra ainsi qu'un seul standard devrait etre 
suffisant pour determiner 1' effectif initial de la 
population d'interet dans 1 1 echantillon d'interet, selon 
un avantage que procure 1' invention. 

Neanmoins, dans une variante qui peut aussi etre prevue si 
necessaire, on obtient les valeurs respectives de mesures 
d'effectifs initiaux F 0st dans plusieurs echantillons 
standard d T effectif s initiaux connus N 0st . On etablit 
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ensuite une dependance entre les effectifs initiaux N 0s t 
des echantillons standard et leurs valeurs respectives de 
mesure d ! effectif initial F 0st . Ensuite, apres la 
determination de la valeur de mesure de l'effectif initial 
pour 1 1 echantillon d'interet F 0 , on determine l'effectif 
initial de la population d'interet N 0 par interpolation 
sur cette dependance. On comprendra que cette dependance 
peut etre encore typiquement une regression du type 
represents sur la figure 5, mais sur laquelle on reportera 
en ordonnees les fluorescences initiales F 0st des standards 
et de 1' echantillon d' interet F 0 (ou leurs valeurs 
logarithmes respectives), plutot que leur valeur d' indic'e 
de basculement C E ep. 

Des lors que 1'on fait appel a un ou plusieurs standards, 
on peut prevoir un ou plusieurs echantillons standard 

d r effectifs initiaux respectifs connus N 0s t, auxquels on 

j 

applique la succession de reactions d' amplifications; 
sensiblement dans les memes conditions que pour 
1 ' Echantillon d'interet. Ensuite, la determination de leur 
valeur de mesure d'effectif initial F 0st se fa'it par 
comparaison des variations parametrees avec les mesures 
exp6rimentales, comme pour 1' echantillon d'interet. 

En d'autres termes, on applique bien entendu les memes 
calculs pour les quantites de fluorescence mesurees et 
ajustees sur le ou les standard (s) et pour 1' echantillon 
d'interet. On estime la quantite de fluorescence F 0sl avant 
le premier cycle, pour le ou les standard (s) selon *le meme 
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procede de determination de F 0 pour 1' echantillon 
d'interet, d6crit ci-avant. 

Un troisieme mode de realisation, correspondant a une 
variante du second mode de realisation ci-avant consiste 
globalement a ajuster le module de l'efficacite E n par 
rapport aux mesures experimentales et a injecter ensuite 
ce modele d' ef f icacite ajuste dans le modele de l'effectif 
courant N n ou encore dans le modele des mesures F n , Ce 
troisieme mode de realisation se resume comme suit. 

La variation parametree construite a 1' etape b) est 
representative du rendement, et, a 1 ' etape c) , on 
determine une variation experimentale du rendement 3 
partir des mesures experimentales, pour comparer la 
variation parametree a la variation experimentale. 
Ensuite, pour acceder a un parametre representa tif de 
l'effectif initial N 0 , a l'etape d) : 

dl) on determine une seconde variation parametree, 
representative de l'effectif courant N n dans 
1 'Echantillon d'interet, et mettant en jeu au moins le 
parametre representant l'indice de basculement C E£P , 
ainsi qu'un parametre. reprfesentatif de l'effectif 
initial N 0 , 

d2) on applique a cette seconde variation une valeur du 
parametre de l'indice de basculement C EEP determine a 
1 1 etape c) , 

d3) et on ajuste au moins le parametre representatif de 
l'effectif initial N 0 par comparaison directe de la 
seconde variation avec les mesures experimentales. 
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On peut avantageusement : 

- a l'6tape d2), appliquer une valeur grossiere de 
l'indice de basculement C E ep, comme on a dfecrit sa 
detection en reference a la figure 4B ci-avant, tandis 
que 

- a l'etape d3) , on peut af finer ensuite la valeur de 
l'indice con j ointement avec 1 1 ajustement du parametre 
representatif de l'effectif initial N 0 . 

Finalement, on retiendra s implement que le second mode de 
realisation des figures 7 et 8, actuellement prefere, se 
distingue de ce troisieme mode de realisation par le fait 
que l'on ne cherche pas a mener les correlations sur : 
1 ' ef f icacite, mais on utilise simplement le modele 
mathematique de variation .de l'eff icacite pour modeliser 
et af finer 1' estimation de la fluorescence corrigee et 
compensee . 

Bien entendu, la presente invention ne se limite pas k la' 
forme de realisation decrite ci-avant a titre d f exemple ; 
elle s'etend a d*autres variantes . 

Ainsi, on comprendra que la presente invention peut 
s' appliquer en outre a la quantification relative, 
notamment par PCR. Dans cette application, on prevoit 
d' amplifier, outre la population d f interet, une population 
de reference soit conjointement dans le meme support, soit 
separement. On relive respectivement : 

- les mesures experimentales representatives de 
l'effectif de la population d ! interet, et 
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les mesures- experimentales representatives de 
l'effectif de la population de reference. 

Le proced6 peut se poursuivre alors par 1 1 application, a 
la population de reference, des etapes a), b) et c) , 
tandis que 1 ' etape d) consiste simplement a determiner un 
ratio entre les effectifs initiaux respectifs de la 
population d'interet et de la population de reference. 

On indique que la quantification relative peut etre 
utilisee pour analyser 1' expression d'un gene d'inter§t au 
cours du developpement d'un organisme. Pour corriger 
notamment les variations de quantite et de quality entre 
les echantillons preleves sur 1' organisme a differents 
instants, on analyse, en plus du gene cible d'interet, un 
gene de reference connu pour avoir un niveau d' expression 
stable au cours du developpement. 

N 

II s'agira finalernent de comparer les ratios 0aM -- entre 

^org- 
ies differents echantillons preleves. 

Pour aboutir aux resultats souhaites, deux strategies sont 
possibles . 



Celle de l'art anterieur basee sur la detection du cycle 
25 Ct se deroule habituellement comme suit. On determine, 
pour chaque echantillon preleve a differents instants tO, 

N 

tl, t2, tn, le ratio 0c ' A/ '' , en utilisant au moins un 

standard (c v est-a-dire un echantillon dont N 0ahl(! et N 0r . f 
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sont connus), ce qui revient a faire deux quantifications 
absolues successives suivies d' un calcui de ratio. 



Une autre strategie particulierement avantageuse dans le 
contexte de l 1 invention consiste a determiner, pour chaque 
echantillon preleve a differents instants tO, tl, t2,>.., tn 



directement le ratio 



^ ^Mf ) ich _ (doable) ech _/Q 



N 



Qablti 



<!ch_lQ 



(^0«/)a/i fO 



\ ^Oref ) 



grace a la formule 



10 Dans cette seconde realisation qui n 1 utilise f inalement. : 
que le parametre Fo, en combinant la technique de' 
1 ! invention, aucun echantillon standard n'est necessaire,f 
de fagon particulierement avantageuse. 



15 On se ref6re maintenant a la figure 10 sur laquelle on a 
representee une installation pour la mise en ceuvre du 
procede selon I 1 invention. Elle comporte un support SUPP 
comportant ici un puits comprenant 1 1 echantillon d'interet 
ECH et un puits comprenant par exemple un echantillon 

20 standard note St. Le support SUPP est enferme dans une 
enceinte ENC par exemple equipee de moyens de chauffage 
(non representes) pour appliquer au standard et a 
1 'echantillon d'interet une reaction de PCR. 
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Dans 1' exemple decrit, on prevoit pref erentiellement 
d'appliquer des mesures de quantities de fluorescence 
emises a chaque cycle par le standard St et par 
l'echantillon d'interet ECH. A cet effet, on introduit 
dans les puits un reactif choisi et on eclaire les 
echantillons a partir d'une lampe (par exemple halogene- 
tungstene) dans le but . de rnesurer des quantit£s 
respectives de fluorescence emanant de l'echantillon 
d'interet et de 1 ' echantillon standard, a chaque cycle de 
PCR qui sera applique. En outre, un appareil de detection 
de fluorescence comporte par exemple un objectif 11 pour 
collecter la lumiere issue de la fluorescence, ainsi que 
des moyens de comptage de photons 10, par exemple une 
camera CCD, et/ou des phqtomultiplicateurs, pour compter 
la fluorescence emise a chaque cycle de PCR, par 
l'echantillon d'interet et le standard. Ainsi, la 
fluorescence emise par chaque puits est avantageusement 
focalisee par 1 ! objectif 11, puis detectee 
preferentiellement par une camera CCD 10 connectee a une 
carte d 1 acquisition 21, par exemple de type PCMCIA prevue 
dans une unite centrale 20 d'un ordinateur. 

Cet ordinateur est alors relie aux moyens de comptage 10 

precites pour recevoir des signaux representatif s de 
mesures de quantite de fluorescence detectee a chaque 

cycle de PCR et traiter ces signaux pour determiner un 

effectif initial de la population d ' interet * avant le 

premier cycle, par la mise en osuvre du procede au sens de 
1' invention. 
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On indique que 1' unite de traitement comporte 
typiquement : 

- une carte d 1 acquisition 21 connectee aux moyens de 
comptage 10, 

- une memoire de travail 25 (par exemple de type RAM) 
pour le stockage temporaire et le traitement des 
signaux pr6cites, 

- une memoire permanente 2 4 pour le stockage du produit 
programme d r ordinateur au sens de l 1 invention et des 
donnees traitees et pretes a exploitation, par exemple 
pour un diagnostic ulterieur, 

- le cas echeant, un lecteur 22 d f un support memoire tel 
qu'un CD-ROM, ou autre, sur lequel peut etre stoc-Jce 
.initialement le produit programme d ' ordinateur , 

- eventuellement, une interface de communication 26 vers 
un site local ou distant (liaison 28), par exemple pour 
la transmission des donnees traitees pour un diagnostic 
sur le patient realise a distance, y. 

- une interface graphique 27 reliee typiquement a : :Un 
moniteur d f affichage 30, 

- et un processeur 23 pour gerer les interactions entre 
ces differents 6quipements. 

On indique que 1* ordinateur peut comporter aussi des 
organes de saisie tels qu'un clavier 41 et/ou une souris 
42 relies a 1' unite centrale 20. 

On retiendra toutefois que 1' installation au sens de 
1' invention comporte globalement : 

- un support SUPP d 1 echantillon au moins pour 
l'echantillon d'interet, 
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- un premier , appareil ENC pour appliquer, au moins a la 
population d'interet dans 1 1 echantillon d'interet, 
ladite succession de reactions d 1 amplification, 

- un second appareil 10 pour relever des mesures 
5 representatives de l'effectif courant de la population 

d'interet, et 

- des moyens inf ormatiques 20 propres a r6cuperer des 
signaux de mesures du second appareil 10 et a mettre en 
oeuvre tout ou partie des.etapes du proc6de au sens de 

10 If invention. 



A cet effet, un produit programme d'ordinateur peut animer 
ces moyens inf ormatiques . Ce programme peut etre stocke 
dans une memoire de 1' unite de traitement 20, ou sur un 

15 support memoire amovible (CD-ROM ou autre) et destine a 
coopSrer avec le lecteur de 1' unite de traitement. Le 
programme d'ordinateur au sens de 1' invention comporte 
alors des instructions pour la mise en ceuvre de tout ou 
partie des etapes du procibde selon 1' invention. Par 

20 exemple, l'algorithme de ce programme peut etre represents 
par un organigramme equivalent au schema de la figure 9. 
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Revendications 

1. Procede mis en ceuvre par des moyens inf ormatiques pour 
quantifier, de facon absolue et/ou relative, une 
population initiale d f acides nucl6iques dans un 
echantillon d f inter§t soumis a une succession 
d 1 applications d ' une reaction' d ' amplification . de la 
population, au cours desquelles on releve des mesures 
experimentales representatives d'un effectif courant de la 
population au moins dans 1 1 echantillon d'int^ret, 

caracterise en ce qu'il comporte les etapes : 

a) prevoir un modele de rendement (E n ) de la reaction 
d ! amplification en fonction de la succession' 
d 1 amplifications, ce modele comportant : 

o une phase sensiblement constante (E 0 ) pour une 
premiere partie des applications de la reaction^ 
et 

o une phase non constante pour une seconde parti§ 
des applications de la reaction, i. 
les premiere et seconde parties etant separees par une 
region . de basculement du rendement entre les phases 
constante et non constante, incluant un indice 
d f application (C EEP ) correspondant sensiblement au 
basculement, 

b) exprimer, en utilisant le modele de rendement, une 
relation faisant intervenir au moins 1' indice de 
basculement et un parametre representant 1' effectif 
initial dans 1 ' echantillon d'interet (N 0 ;F 0 ), 

c) par comparaison avec les mesures experimentales, 
determiner au moins l f indice de basculement, et, dans 
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une etape d) ulterieure ou conjointe, en deduire 
l'effectif initial dans 1 1 echantillon d'interet (N 0 ) . 

2. Procede selon la revendication 1, caracterise en ce 
que : 

- a l'6tape b) , on exprime, en utilisant le modele de 
rendement, une varialion parametree en fonction de 
ladite succession d 1 amplifications, mettant en jeu au 
moins un parametre representant l'indice de basculement 
(Ceep) f et « 

- a l'6tape c) , par comparaison avec lesdites mesures 
. experimentales, on determine au moins ce parametre 

representant l'indice de basculement (C EBP ) . 

3. Procede selon la revendication 2, caracterise en ce que 
ladite variation parametree est representative de 
l'effectif courant (N n ) dans 1 1 echantillon d'intSret, 

en ce que- ladite variation met en jeu en outre un 
parametre representatif de l'effectif initial dans 
1 1 echantillon d'interet (N 0 ) , 

et en ce que, aux etapes c) et d) , on determine 
sensiblement con jointement les parametres representatif s 
dudit indice d 1 application (C EEP ) et de l'effectif initial 
(N 0 ) dans 1 ' echantillon d'interet. 

4. Procede selon la revendication 2, caracterise en ce que 
ladite variation parametree est representative du 
rendement, et en ce que, a 1' etape c) , on determine une 
variation experimentale du rendement a partir desdites 
mesures experimentales, pour comparer la variation 
parametree a la variation experimentale. 
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5. Procede selon la revendication 4, caracterise en ce 
que, a l'etape d) : 

dl) on determine une seconde variation parametree, 
representative de l'effectif courant (N n ) dans 
1 ' echantillon d'interet, et mettant en jeu au moins le 
• parametre- representant ledit indice d 'application 
(Ceep) , ainsi qu'un parametre representatif de 
l'effectif initial (N 0 ) dans 1 ' echantillon d'interet, 

d2) on applique a la seconde variation une valeur du 
parametre de 1' indice (C EEP ) determine a l'etape c),. 

d3) et on ajuste au moins le parametre representatif de 
l'effectif initial " (N 0 ) par comparaison directe de £la 
seconde variation avec les mesures experimentales. 

6. Procede selon la revendication 5, caracterise en ,,ce 
que : 

£ l'etape d2), on applique une valeur grossi&re -,de 
1' indice (C EEP ) , tandis que '4 

- a l'etape d3), on affine la valeur de l'indice 
conjointement avec l'ajustement du parametre 
representatif de l'effectif initial (N 0 ) . 

7. Procede selon la revendication 3 ou la revendication 6, 
caracterise en ce que ladite variation parametree ou, 
respectivement, ladite seconde variation parametree : 

- est representative desdites mesures experimentales, et 
— comporte un parametre cor respondant a une valeur de 

mesure (F 0 ) representative de l'effectif initial, 
et en ce que 1 ' on determine la valeur de mesure de 
l'effectif initial (F 0 ) par comparaison desdites 
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variations parametrees (F n ) 
experimentales . 



avec 



les 



mesures 



8. Procede selon l ,f une des revendications precedentes, 
caracterise en ce qu f on applique une etape prealable de 
traitement des mesures experimentales comprenant une 
soustraction d'un * bruit de fond - de mesures et 
d' introduction d'une compensation pour tenir compte d'une 
mesure non nulle (e) representative de 1' effectif initial. 

9. Procede selon l ! une des revendications 7 et 8, 
caracterise en ce qu 1 on obtient une valeur de mesure d'un 
effectif initial (F 0 3t ) dans un echanzillon standard 
d'effectif initial connu (N 0 st ) , 

et en ce qu ' on determine la valeur de 1» effectif initial 
dans l'6chantillon d'interet (N 0 ) par application d'une 
meme relation de proportionnalite entre 1' effectif initial 
et sa mesure, dans 1 ' echantillon standard d'une part, et 
dans 1' echantillon d'interet d'autre part. 

10. Procede selon 1 ' une des revendications 7 et 8, 
caracterise en ce qu 1 on obtient des valeurs respectives de 
mesures d r effectifs initiaux (F 0 st ) dans des echantillons 
standards d'effectifs initiaux connus (N 0 st ) , 

et en ce que : 

- on etablit une dependance entre les effectifs initiaux 
(N 0 3t ) des echantillons standard et leurs valeurs 
respectives de mesure d'effectif initial (F 0 3t ), 

- et, apres la determination de la valeur de mesure de 
l 1 effectif initial pour 1 1 echantillon d'interet (F 0 ) , on 
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determine l'effectif initial de la population d'interet 
(N 0 ) par interpolation sur ladite dependance. 

11. Procede selon 1 ' une des revendications 9 et 10, 
caracterise en ce que 1 1 on prevoit un ou une pluralite 
d' echantillons standard d'effectifs initiaux respectifs 
connus (N 0 st ) , auxquels on applique ladite succession 
d ' amplifications, sensiblement dans les memes conditions 
que pour 1 ' echantillon d'interet, ainsi que la 
determination de leur valeur de mesure d'effectif initial 
(F 0 st ) par comparaison des variations parametrees (F n 5t ) 
avec les mesures experimentales . 

12. Procede selon la revendication 3, caracterise en.rce 
que l f on prevoit une pluralite d T echantillons standard 
d ! effectifs initiaux respectifs connus, auxquels .> on 
applique ladite succession d ! amplifications, sensiblement 
dans les memes conditions que pour 1 ' echantillon 
d'interet, ainsi que, la determination de leurs indices 
respectifs selon les etapes a), b) et c) , 

et en ce que, a l'etape d) : 

- on etablit une dependance entre les effectifs initiaux 
des echantillons standard (N 0 st ) et leurs indices 

(C E £p St ) / 

- et, apres la determination de l'indice pour 
1 1 Echantillon d'interet (C E ep) / on determine l'effectif 
initial de la population d'interet (N 0 ) par 
interpolation sur ladite dependance. 

13. Proced6 selon l'une des revendications 1 a 8, 
caracterise en ce que, pour une quantification relative, 
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on pr^voit, outre la population d'interet, une population 
de reference soumise a une succession d 1 applications de 
reaction d 1 amplification, et on releve respectivement : 

- les mesures experimentales representatives de 
l'effectif de la population d'interet, et 

- les mesures experimentales representatives de 
l'effectif de la population de reference, 

le procede se poursuivant par 1 ' application, a la 
population de reference, des etapes a), b) et c) , tandis 
que l'etape d) consiste a determiner un ratio entre les 
effectifs initiaux respectifs de la population d'interet 
et de la population de reference. 

14. Procede selon l'une des revendications precedentes, 
cafacterise en ce que : 

- on exprime les mesures experimentales sous la forme 
d'une variation experimentale du rendement en fonction 
de ladite succession d ' applications , 

- on obtient une variation experimentale du rendement en 
fonction de ladite succession duplications, 
comportant : 

* une premiere region sensiblement bruitee pour des 
indices d ' applications faibles, 

* suivie d'une seconde region moins bruitee pour des 
indices d ' applications plus eleves. 

15. Procede selon la revendication 14, dans lequel ladite 
phase non constante du rendement est decroissante, 
caracterise en ce que : 

- on estime une valeur grossiere de phase constante du 
rendement (E 0 ) , et 
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- au .moins pour la recherche de l 1 indice (C E ep) dans 
ladite r6gion de basculement, on ignore une partie au 
moins des mesures dans ladite seconde region moins 
bruitee et dont le rendement estime est inf6rieur a une 
valeur seuil, de preference a une fraction de la phase 
constante (Eo) . 

16. Procede selon 1 1 une des revendications 14 et 15, dans 
lequel ladite phase non constante du rendement est 
decroissante, caracterise en ce que 1 ' on rep£re ladite 
region de basculement, dans le sens des indices 
d' applications decroissants a partir de ladite seconde 
region moins bruitee, en detectant un indice grossier 
d' application d 1 amplification dont le rendement depasse 
sensiblement une valeur predetermine (E 0 =l) . 

17. Procede selon la revendication 16, caracterise en ce 
qu'on affine l 1 estimation de la valeur, le cas ech:eant 
fractionnaire, . dudit indice d 1 application (C EEP ) dan$' la 
region de basculement, dans le sens des indices 
d ' applications croissants a partir de 1 1 indice grossier, 
en detectant un indice d r application d ' amplification dont 
le rendement vaut environ ladite valeur predeterminee 
(E 0 =l) • 

18. Procede selon l'une des revendications precedentes, 
caracterise en ce que ladite phase -non constante du 
rendement est modelisee par une exponentielle decroissante 
comportant un parametre de decroissance (p), et en ce que 
ledit parametre de decroissance (P) est determine, a 
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l'etape c) avec l ! indice dans la region de basculement 
(Cbep) / par comparaison avec les mesures experimentales . 

19. Procdde selon 1 ■ une des revendications precedentes, 
5 dans lequel la reaction d 1 amplification est une reaction 

de PCR ("Polymerase Chain Reaction") en temps reel. 

20. Procede selon l'une des revendications precedentes, 
dans lequel lesdites mesures sont des quantites de 

10 fluorescence emises. 



15 



21. Installation pour la mise en ceuvre du procede selon 
l'une des revendications precedentes, caracterisee en ce 
qu'elle comporte : 

- un support d r echantillon au moins pour 1 1 echantillon 
d'intferet, 

- un premier appareil pour appliquer, au moins a la 
population d ' interet dans 1 ' echantillon d ' interet , 
ladite succession de reactions d 1 amplification, 

20 " un second appareil pour relever, des mesures 
representatives de l'effectif courant de la population 
d'interet, et 

- des moyens inf ormatiques propres a recuperer des 
signaux de mesures du second appareil et a mettre en 

25 oeuvre tout ou partie des etapes du procede selon l'une 

des revendications precedentes. 

22. Produit programme d * ordinateur , destine a etre stocke 
dans une memoire d'une unite.de traitement, ou sur un 

30 support memoire amovible et destine a coop6rer avec un 
lecteur de ladite unite de traitement, caracterise en ce 
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qu'il comporte des instructions pour la mise en oeuvre de 
tout ou partie des etapes du proc^de selon l'une des 
revendications 1 a 20. 
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